Metallic behavior in Si/SiGe 2D electron systems 
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We calculate the temperature, density, and parallel magnetic field dependence of low temperature 
electronic resistivity in 2D high-mobility Si/SiGe quantum structures, assuming the conductivity 
limiting mechanism to be carrier scattering by screened random charged Coulombic impurity cen- 
ters. We obtain comprehensive agreement with existing experimental transport data, compellingly 
establishing that the observed 2D metallic behavior in low-density Si/SiGe systems arises from the 
peculiar nature of 2D screening of long-range impurity disorder. In particular, our theory correctly 
predicts the experimentally observed metallic temperature dependence of 2D resistivity in the fully 
spin-polarized system. 
PACS Number : 71.30.-f-h; 73.40. Qv 
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The two-dimensional (2D) "metallicity" is by now a 
well-established ubiquitous low-temperature phenomena 
[1] where a relatively high-quahty (i.e. low disorder) 2D 
carrier (either electron or hole) system confined in semi- 
conductor heterostructures exhibits a pronounced tem- 
perature (T) and/or parallel (i.e. in the 2D plane) mag- 
netic field (S||) dependent resistivity p{T; i?||) in the di- 
lute or low carrier density regime. The strength of the 
metallicity, i.e. the relative magnitude of the variation 
of p{T) with T, is found to depend on the specific 2D 
material involved as also does the precise low-density 
regime where the strong metallicity manifests itself. For 
example, in 2D Si MOSFETs the metallicity occurs [2] 



around 10^ 



with p{T) changing by as much as 



a factor of 3 — 4 as temperature changes from 50 mK 
to 5K (at B = 0) whereas in 2D n-GaAs electron sys- 
tems the metallicity is observed [3] for density well be- 
low lO^^CTO^^ with p{T) changing only by about 25% as 
T changes from 50 mK to a few K. The 2D metallicity 
has been experimentally observed in essentially all exist- 
ing semiconductor-based 2D systems, most notable (and 
most extensively studied) being electrons in Si MOS- 
FETs, holes in p-GaAs 2D structures [4] and very re- 
cently, electrons in n-GaAs 2D structures [3]. 

One of the more important systems to manifest 2D 
metallic behavior is the 2D electron system in modula- 
tion doped Si/SiGe structures [5-8] where the electrons 
are confined by the Si-SiGe interface potential barrier 
on the Si side of the heterostructure and the modula- 
tion doping delta layer (contributing the 2D electrons) 
is typically placed 100 - 200 A inside the insulating 
SiGe side. The relatively impurity-free high quality of 
the materials (both Si and SiGe), and the lack of in- 
terface impurities and roughness (at the Si-SiGe inter- 
face) make the Si/SiGe 2D electron system essentially 
ideal from the perspective of understanding 2D metallic 
behavior. In particular, the dominant low-temperature 
(below 5K or so where phonon scattering should be ex- 
ponentially suppressed in silicon) scattering mechanism 
in Si/SiGe 2D electron system is thought to be scatter- 
ing by the remote dopants in the modulation delta layer 



which, being far (i.e. 100 - 200 A) away from the 2D elec- 
trons residing inside Si, leads to very high, by more than 
an order of magnitude higher than the corresponding Si 
MOSFETs, carrier mobilities. The high quality of the 
system and the simple isotropic parabolic 2D energy dis- 
persion (in contrast to the p-GaAs based 2D hole systems 
which suffer from the multiband valence band structure, 
strong 2D anisotropy, strong spin-orbit coupling, inter- 
spin-split-band scattering, and strong phonon scattering 
leading to considerable complications in its 2D metallic- 
ity) make Si/SiGe based 2D electron systems an excellent 
candidate for studying 2D metallic behavior. In particu- 
lar, 2D electrons in Si/SiGe structures combine, from the 
2D metallic physics perspective, the advantages of n-Si 
MOSFETs (i.e. strong metallicity) with the advantages 
of the n-GaAs system (i.e., high mobility and low den- 
sity) without the n-GaAs disadvantage of weak metal- 
licity and phonon scattering problems. Metallicity has 
been observed [8] in 2D Si/SiGe systems down to densi- 
ties as low as 3 X lO-'^^cm"^, and the relevant interaction 
parameter (r^) is typically 5-10 for this system. 

In the context of 2D metallic physics it is therefore 
of considerable interest that recently two independent 
groups have reported detailed experimental data on the 
temperature, density, and parallel magnetic field depen- 
dence of 2D resistivity in Si/SiGe confined electron sys- 
tems. The reported Si/SiGe experimental results, partic- 
ularly in the presence of a parallel magnetic field, have 
significant qualitative difference with the corresponding 
experimental results in n-Si MOSFETs, warranting a the- 
oretical investigation. We provide, in this Communica- 
tion, such a transport theory for 2D electrons in Si/SiGe 
structures, finding good agreement between our theoret- 
ical and the existing experimental results. 

We use a minimal zeroth-order theoretical model [9] 
to describe the 2D transport properties, keeping the to- 
tal number of unknown parameters to a minimum. We 
neglect all phonon scattering effects although it is fairly 
straightforward to include them, mainly because our the- 
oretical estimate shows phonon scattering to be negligi- 
ble for 2D electrons in Si/SiGe structures in the T < 5K 
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regime of interest to us. We assume that the 2D car- 
rier conductivity is entirely limited by screened impu- 
rity scattering, where the disorder arises from random 
background charged impurity centers (i.e. unintentional 
dopants) and the ionized dopants (randomly distributed) 
in the modulation dopant layer. The background unin- 
tentional impurity density is extremely low (^ 10^^cm~^ , 
which is equivalent to a 2D density of 10*cm~^), as im- 
plied by the high achieved mobilities (~ 7 x lO^cm^ /Vs) 
in gated Si/SiGe quantum well structures. The 2D resis- 
tivity is given in the Boltzmann theory by p = m/ne^ir), 
where n,m,T are respectively the 2D carrier density, the 
carrier effective mass, the transport relaxation time, and 




with /(e) being the Fermi distribution function and 

X (1 - cos6'kk')'5(ek - ek')- (2) 

Eq. (1) indicates a thermal average over the energy de- 
pendent relaxation time T(e), with e = k"^ /2m being 
the usual parabolic 2D electron energy dispersion. In 
Eq. (2), u^~^{q;z) is the (finite-temperature) screened 
electron-impurity interaction for 2D momentum trans- 
fer of q and Ni{z) is the random charged impurity den- 
sity in the direction [z) normal to the 2D {x-y) plane 
of confinement of the electron layer (the 2 = being 
the Si/SiGe interface plane), which is given in our model 
by Ni{z) = + ni6{z - z,), where Nf^ is the back- 
ground 3D charged impurity density and is the 2D 
charged dopant density in the modulation doping (5-layer 
at z = Zi. We take Zi to coincide with the modulation 
doping layer in the experimental Si/SiGe samples. The 
key quantity in Eqs. (l)-(2) is the screened charged im- 
purity disorder potential which we take to be 

u'-\q,z)=v'-\q;z)/s{q), (3) 

where u*^"* is the bare Coulomb potential for electron- 
charged impurity interaction, and e{q) = 1 — v{q)Il{q) 
is the RPA dielectric screening function due to the 2D 
electrons themselves, where v{q) is the 2D bare Coulomb 
electron-electron interaction, and Il(q) = Il{q; T) is the 
2D finite-temperature and finite wave vector polarizabil- 
ity function. All our calculations include the realistic 
quasi-2D quantum form factor effects, arising from the fi- 
nite width of the 2D layer in 2-direction, which enter the 
expressions for v{q) and i/^'((jr), that are of considerable 
quantitative importance. All quantities entering the the- 
ory are known except for Nf^ and n^, which are the only 
two parameters of the model. We set Ni = 2x 10^'^cm~^, 
rii = 3.5 X 10^°cm~^, Zi = lOOA throughout our calcu- 
lations, instead of using them as adjustable parameters 



since onr interest here is a qiialitative (and perhaps semi- 
quantitative) understanding of 2D metallicity in Si/SiGe 
systems. Note that there is little scattering by inter- 
face impurities and interface roughness in high-mobility 
Si/SiGe 2D systems, which is their distinguishing feature 
compared with Si MOS systems, leading to the extraor- 
dinarily (a factor of 20 - 40 higher than in Si MOSFETs) 
high 2D mobility. 
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FIG. 1. The calculated 2D mobility of Si/SiGe 2D systems 
at 0.3K as a function of carrier density for different values of 
Zi. Dashed line represents ji oc n°'®. 

In Fig. 1 we show our calculated 2D mobility, /x = 

1/nep, as a function of 2D electron density at T = 0.3K 
for a number of different values of Zi, the position of the 
2D dopant layer. The calculated mobility for Zi = 100 A 
agrees well (quantitatively for n > 0.5 x IQ^^crn,^^) with 
the density-dependent mobility measured experimentally 
[8] which also (in agreement with Fig. 1) manifests a 
H oc nP'^ dependence on density at higher densities in 
contrast to /it oc behavior observed in Si MOSFETs 
(i.e., our Zi = result corresponding to scattering by 
interface impurities) . This establishes the validity of our 
model for the experimental Si/SiGe 2D systems. 

In Fig. 2 - 4 we show our calculated results for the 
2D resistivity in Si/SiCc systems. In Fig. 2 wc show 
the zero-field (i^y = 0) temperature dependence of the 
2D resistivity for a number of carrier densities in the 
2D metallic phase. The theoretical results are in ex- 
cellent qualitative agreement with experimental obser- 
vations [5-8] with the maximum temperature induced 
fractional change in the resistivity, Ap/po where p = 
po + ApiT) with pq = p{T — > 0), being about 2 at 
lower densities in the T « lOOmK — 5K range. At 
higher densities (> lO-'^-'^cm"^), Ap/Ag is, consistent with 
the experimental data, only about 50% (in the — 5K 
range). The relatively weaker temperature dependence 
in Si/SiGc 2D system compared with Si MOS system 
arises entirely from the suppressed role of 2fci?-scattering 
in the modulation doped Si/SiGe systems - in fact, the 
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Si/SiGe system will manifest very strong (stronger than 
Si MOSFET) 2D metallicity if its mobility would be lim- 
ited by interface scattering as in Si MOSFETs (or purely 
by the background impurities as in p-GaAs and n-GaAs 
2D systems). It is well-known that the strong metallicity 
in the transport properties of low-disorder 2D systems 
arises from the temperature induced suppression of the 
2fcF-Kohn anomaly in screening (i.e. the T = cusp at 
q = 2kF in the 2D polarizability function), and the over- 
all strong suppression of 2kF scattering in the modulation 
doped structures [due to the presence of the e"^'^^^'!^*! fac- 
tor in the electron- impurity bare interaction v^~'^{q; z) in 
Eq. (3)] leads to relatively weaker 2D metallicity. In 
fact, the temperature dependence would be substantially 
weaker than that in Fig. 2 (and therefore substantially 
weaker than the experimental observation) if we set the 
background impurity density Nf^ to be zero in our calcu- 
lation. The temperature dependence (Fig. 2) of p{T; n) 
is thus a simple diagnostic for the strength of background 
impurity scattering in Si/SiGe system just as the density 
dependence (Fig. 1) is a diagnostic for the remote im- 
purity scattering. We also mention the "almost parallel" 
behavior of p{T) at various densities in Fig. 2, which has 
been experimentally seen [8] . 
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FIG. 2. Calculated resistivity as a function of temperature, 

p(r), of 2D Si/SiGe systems for various densities n =4.2, 5.15, 
6.3, 8.2, 10.8 14.5xl0^°cm"^ (top to bottom). 

In Fig. 2 p{T) shows a nonmonotonicity (particularly 
prominent at lower densities) where it increases with T 
at first and then decreases with increasing T after going 
through a shallow maximum at a characteristic density- 
dependent temperature Tp{n) where dp/dT changes sign, 
i.e. dp/dT = at T = Tp. Our theoretical Tp{n) in Fig. 2 
scales approximately with the Fermi temperature Tp , but 
is not identical to Tp. This is precisely the experimental 
observation [8]. This non-monotonicity, which is obvi- 
ously a non- asymptotic finite temperature phenomenon, 
arises from the competition between an increasing and 
decreasing contribution to p{T) as defined by Eqs. (1) 



- (3). The temperature dependence of the dielectric 
screening function in Eq. (3), with screening (partic- 
ularly the strongly resistive 2fci?-scattering) decreasing 
with increasing temperature, gives rise to increasing ef- 
fective disorder with increasing T, and hence increasing 
p{T) with T. At higher temperatures the thermal aver- 
aging in Eq. (1), which always leads to p{T) decreasing 
with T, becomes quantitatively more important, giving 
rise eventually to a p{T) decreasing with T for T > Tp. 
In the classical non-degenerate regime the thermal av- 
eraging always wins out and impurity scattering limited 
p{T) necessarily decreases with increasing T (which is 
why we earlier called this nonmonotonicity somewhat 
simplistically "the quantum-classical crossover", gener- 
ating some confusion on this matter), but thermal aver- 
aging always (i.e. even in the quantum T <C Tp regime) 
leads to a negative sub-leading contribution to p{T) going 
as 0{T/Tp)^, whereas the screening contributes the lead- 
ing order term going as 0{T/Tp), so in the T/Tp 
limit the increasing contribution arising from 2k p screen- 
ing wins out. Depending on the details of the system, 
Tp[n) could be quite low (e.g. much lower than Tp in 
weakly screening systems) although it should scale ap- 
proximately with Tp. Our theoretical Tp{n) in Fig. 2 is 
roughly of the order of the experimental Tp, but the the- 
oretical Tp is somewhat higher quantitatively than the 
experimental Tp. Many details (e.g. phonon scattering, 
exact charged impurity distribution, the confinement po- 
tential, the effective mass, etc.) affect the precise value 
of Tp, and therefore a quantitative comparison between 
theory and experiment is not particularly meaningful. 

In Figs. 3 and 4 we show our calculated results for 
the 2D magnetoresistance in Si/SiGc system in the pres- 
ence of a parallel magnetic field which polarizes elec- 
tron spins, consequently weakening screening consider- 
ably (i.e. Bji acts similar to finite T effects). Of the three 
physical mechanisms [10] that affect 2D parallel field 
magnetoresistance, we have ignored the magneto-orbit 
effect [11] because it is negligibly small in the Si system 
at the applied parallel field values of experimental inter- 
est. We include the other two (competing) parallel field 
corrections, both arising from the electronic spin polar- 
ization induced by the applied field (and therefore both 
saturating for _B|| > Be where the system is completely 
spin polarized) these being the spin-polarization in- 
duced weakening of screening (leading to positive magne- 
toresistance) and the spin-polarization induced enhance- 
ment of the 2D Fermi wave vector kp (leading to negative 
magnetoresistance). The Si 2D system of interest to us 
is a strongly screening system with qTp/2kp ^ 1, where 
qpp is the 2D screening wave vector, and therefore the 
screening effect is by far the dominant parallel field effect, 
producing a strong positive magnetoresistance apparent 
in Figs. 3 and 4 (and in the experiments). 
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FIG. 3. Calculated magnetoresistance as a function of 
in-plane magnetic field for different values of the impurity sep- 
aration Zi — 0, 50, 100 A at a fixed 2D density n = 10^^ cm~^ 
and at T = 0.3K (solid lines) and 3. OK (dashed lines). 

In Fig. 3 we show the calculated magnetoresistance 
/?(iJ||)//9(0) for different values of the impurity separa- 
tion Zi. As expected, larger values of zi drastically reduce 
the maximum magnetoresistance since the separation be- 
tween the impurities and the 2D electrons suppresses the 
quantitative importance of 2kF resistive scattering. Our 
maximum theoretical value of p{B^^) / p{0) for Zj = lOOA 
of about 1.7 is in excellent agreement with the corre- 
sponding experimentally observed [8] value, lending sup- 
port to our theoretical screening picture. 

In Fig. 4, which is perhaps the most important result 
presented in this work, we show our calculated magne- 
toresistance as a function of (at a fixed 2D density of 
lO^^cm^^) for several different temperatures - the inset 
of Fig. 4 shows p{T) at B\\ = and 9T (> Be ~ 5T 
at lO^^cmr^ density). The theoretical results in Fig. 4 
are qualitatively similar to the corresponding experimen- 
tal results. The most important feature of Fig. 4 is the 
metallic temperature dependence oi p{T; _B||) at i?|| > Be, 
which, although suppressed from the strong zero-field 
metallic temperature dependence p(T; B|| = 0) due to 
the spin-polarization induced suppression of screening as 
the spin degeneracy changes from the paramagnetic value 
of 2 for B|| = to the "ferromagnetic" value of 1 for 
^11 > -Be, is still strongly metallic in character with p{T) 
changing by almost 50 % at B\\ = 9T « 2Bc. This is 
precisely the experimental observation. In fact, our cal- 
culated temperature dependence of p{T; _B|| > Be) is in 
excellent qualitative agreement with the experimental re- 
sults [7]. 

The significance of our theoretical agreement with the 
observed "metallic" (i.e. dp/dT > 0) resistivity for 
> Be arises from the very specific and contrasting 
prediction for p{T: _B|| > Be) in the so-called "interaction 
theory" [12], which purports to extend the screening the- 
ory by calculating the electron-electron interaction effects 



(to leading order in T/Tp) on the 2D transport prop- 
erties to all order (perturbatively) in interaction for an 
unrealistic zero-range white-noise model of bare impurity 
disorder. The categorical prediction [12] of the interac- 
tion theory is that for a fully spin polarized 2D system 
(i.e. for B|| > Be) the 2D resistivity must necessarily 
manifest an "insulating" (i.e., dp/dT < for _B|| > Be) 
behavior. This interaction theory prediction obviously 
disagrees qualitatively with our results in Fig. 4, where 
dp/dT > for B\\ > Be, as well as with experimental 
observations [7] in the Si/SiGe based high- mobility 2D 
electron systems. This qualitative disagreement (agree- 
ment) between interaction (screening) theory and exper- 
iment shows that the interaction theory is invalid in the 
2D Si/SiGe systems due to the long-range nature of the 
underlying impurity disorder. We have also calculated 
(not shown) p(i?||) as a function of at various densi- 
ties, finding that p(B||)/p(0) at various densities approx- 
imately scale with B\\/Be as reported experimentally [8]. 
But we find that this scaling is non-universal and some- 
what temperature dependent. 
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FIG. 4. Calculated magnetoresistance as a function of 
in-planc magnetic field for different temperatures at a fixed 
carrier density n = lO'^'^ cm~^. Inset shows the resistivities 
p(T, B\\ ) at B\\ = and 9T as a function of temperature. 

In summary, we have developed a realistic screening 
theory for the density, temperature, and parallel mag- 
netic field dependence of low-temperature 2D resistivity 
p(r;n;i?||) in Si/SiGe electron systems. Our results are 
in comprehensive qualitative agreement (in T, n, and By 
dependence) with recent experimental observations [5 8]. 
In particular, we get agreement with the experimental 
observation of a "metallic" (i.e., dp/dT > 0) 2D resistiv- 
ity in the completely spin-polarizcd (i?|| > Be) system, 
which disagrees qualitatively with the opposite prediction 
(i.e. dp/dT < for B|| > Be) of the interaction theory 
[12]. The agreement between our theory and experiments 
directly demonstrates the key role of long-range Coulom- 
bic disorder in 2D metallicity. 
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